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For a broad class of unitary ensembles of random matrices we demonstrate the universal 
nature of the Janossy densities of eigenvalues near the spectral edge, providing a different for- 
mulation of the probability distributions of the limiting second, third, etc. largest eigenvalues of 
the ensembles in question. The approach is based on a representation of the Janossy densities 
in terms of a system of orthogonal polynomials, plus the steepest descent method of Deift and 
Zhou for the asymptotic analysis of the associated Riemann-Hilbert problem. 



1 Introduction 

. Consider the probability measure P n on the space ofnxn Hermitian matrices M defined by 



dP n {M) = _L e -«trV(M) ^ 
An 



in which tr denotes the matrix trace, dM is the Lebesgue measure, and the potential V grows 
sufficiently fast at ±oo so that the normalizer Z n < oo. This prescription is an instance of the 
unitary ensembles of Random Matrix Theory; the invariance dP n {U* MU) = dP n (M) for any nx n 
unitary matrix U explains the terminology. 

Regarding their spectral properties these ensembles are integrable. That is to say, the joint 
probability density of the eigenvalues x\, X2, ■ ■ ■ , x n induced by P n may be computed: 

p n (x 1 ,...,x n ) = ^ r J] \x e -x k \ 2 e- n ^ v ^\ (1.1) 

l<£<k<n 

with a new normalizer Z n . Even more, all finite dimensional correlation functions of the eigenvalues, 

I roc r-oo 

p^\x l ,...,x k ) = — / •••/ p n (xi,...,x k ,x k+ i,...,x n )dx k+1 ---dx n , (1.2) 

have explicit expressions. Bring in the system of polynomials 

Pk,n( X ) = lk,nX k + . . . , 
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k = 1, . . . , n with 7fc >n > 0, orthonormal with respect to the weight w n (x) = e nV ( x ) over R. That 
is, JT 00 Pe,n(x)Pk,n(x)wn(x)dx = and it holds 

p^(xi, ■ ■ ■ ,x k ) = det K n (xi,x m ) , (1.3) 

J l<£,m<k 

in which 

n-l 



K n {x,y) = \Jw n {x)\Jw n {y) } J Pk,n(x)pk, n (y) (1-4) 

fc=0 

/ / rr7n-l,n Pn,n(^)Pn-l,n(y) -Pn-l,n(^)Pn,n(y) 



7n,n X - V 



by the formula of Christoffel-Darboux. The form of (jl.3p implies that the ensemble eigenvalues 
comprise a determinantal point process. 

With the above normalization, pn (xi, . . . , Xk) is really a joint intensity of there being an 
eigenvalue, irrespective of order, at each of the points x\ through x^. Alternatively, fix a subset V 
of M containing xi, . . . ,x^. Then, the probability that there are exactly k eigenvalues in V, one at 
each of those same points, defines the fc-th level Janossy density, denoted by J^{x\, . . . ,Xk)- For 
any determinantal point processes the Janossy densities are also determinantal ([5] p. 140): in our 
case, 

jj, k }(xi, . . . ,x k ) = D(T) x det L n r(xe,x m ) , (1.5) 

L J l<£,m<k 

where 

L n>r = K n>T (I-K n>r )- 1 , (1.6) 

the kernel K n ^(x,y) equaling lr(x)K n (x,y)lr(y), and the prefactor D(T) is the Fredholm deter- 
minant 

D(T) = det(I - K n! r). (1.7) 

More important for what follows, it has recently been shown in [2] that kernel of L n p is also 
Christoffel-Darboux type. In particular, 

T , \ I T\ I mV 1 .™ Pn,n( X )Pn-l,n(y) ~ Pn-l,n{ X )Pn,n{v) h „^ 

L n ,r{x,y) = y/w n (x)y/w n {y) — — ■ '■ '■ '■ '■ , (1.8) 

7n,n X-y 

where {pk,n} are the polynomials orthogonal to the weight w n (x), now restricted to the complement 
of T: 

Pe,n(x)Pk,n(x)w n (x)dx = 5 ik - 



\r 

For a large class of potentials V, |11| employs the Riemann-Hilbert Problem (RHP) characteriza- 
tion of the system {p k ,n} to obtain sharp n— >oo asymptotics of the kernel K n , and thus also the 

(k) 

correlation functions p n , in the bulk of the spectrum. Here we take up the analogous project for 
the Janossy densities at the spectral edge by analyzing the {pk,n} system. 
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The requirements on the potential V are described in terms of the equilibrium measure fiy, or 
weak limit of the eigenvalue counting measure. This may be characterized as the infimum of 

/OO yOO /"OO 
/ ^g-——rdfx(x)dfi(y)+ V(x)d^x), (1.9) 
-oo J —oo |*" y\ J — oo 

over the space of probability measures on R. Now, if 

V : R -> R is real analytic (1.10) 

and 

lim V ^ - = +00, (1.11) 

IxHoo log(a; 2 + 1) v ; 

then [7] proves that this infimum is uniquely attained at fly. Further, \iy possesses a density ipy(x) 

with compact support comprised of a finite number of intervals. Assumptions fjl . lOf) and (II. lip are 

adopted here. By a scaling we may fix the rightmost edge of the support of ipy(x) at x = 1, and 

we further assume that 

iftv( x ) is regular. (1-12) 

By this we will mean the following. 

(a) Vv vanishes like a square-root at each endpoint of supp(/uy). 

(b) Vv is strictly positive in the interior of supply)- 

(c) Strict inequality holds in the characterizing Euler-Lagrange equations in the 
exterior of supp(^v)> see ()2.10p . 

One imagines that square-root vanishing at x = 1 would suffice; full regularity has been assumed 
for technical reasons. 

For V(x) satisfying (jl.lUp . (|1.11|) and (|1.12p . one may infer from the results in [11] that the 
kernel Kjy at the spectral edge has the universal limit, 

1 v , x n , y \ _ Ai(x)Ai'(y) - Ai'(x)Ai(y) 



lim m K n [l + ^,1 + ^7^ )=^ J — (1-13) 

with constant cy > and Ai(-) the Airy function^ Based on this, it is expected that the kernel 
L W) r for r = [1 + a/(cyn 2 / 3 ), 00) with any real a will have a universal limit as re^oo. We introduce 
the shorthand, 

L n ,a(x,y) = L n,ll+ a /(c v nV3 )tO0) {x,y), (1.14) 

and, noting that the regime of interest is for x, y 6 V, prove the following. 

Theorem 1.1. Assume that the potential V{x) satisfies il.lO\) . U.ll\) and Then, there 

are pairs of functions {fa( z )>9a( z )} an d {fa( z )i9a( z )} defined for a > and a < respectively, 
such that the following universal asymptotics hold. For a > 0, 

1 r L w , a (i+ * 1+* ) = f^iy)-9Mf-{y) +0(n - 2/3)> (L15) 



cyn 2 / 3 ' V cyn 2 / 3 ' cyn 2 / 3 / x — y 



1 Though understood to hold in greater generality, a detailed proof of (|1.13|) actually only appears in the literature 
for polynomial V [10] , 
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while for a < 0, 

1 T (-, , x t , V \ _ fa( x ~ a )9g(y -a)~ gg(x - a)fg(y - a) , _ 2 /3n 
cyn z > 6 V cyn 1 ' 6 cyn 1 ' 6 ) x — y 

(1.16) 

-Bot/i estimates are uniform for x and y restricted to compact sets of (a, oo). 

The functions f&(z) and g^{z) are read off from the solutions of a pair of 2 x 2 RHPs denoted 
by RHP^ for a > and RHP^ for a < 0. With the corresponding contours and their orientations 
depicted in Figure [Q, we have: 

RHP^ (a > 0): Seek a 2 x 2 matrix valued function M~*(z), analytic in C\£~* such that: 

_4 2 3/2 \ 

(.U-) +l :)-(.U-)_(:)j ^ 6 ^ J, z€(0,a), 

+ - (-V-)-C) ( 4 J 3/2 M, argz = ±§7r, (1.17) 
e3 z 1 J 

(M-) + (z) = (M-) + (z) z ^ (-°°'°) 

with 

AT*(z) =z~\ a3 f | | j e-x CT 3(j + C )( z - 1 )) ) z ^oo. (1.18) 
RHP'" (a < 0): Now seek a 2 x 2 matrix valued function M~*(z), analytic in C\S^ such that: 



(AT)+(z) = (AT)_(z) I e | z 3/2 +2az i/ 2 1 j, argz = ±§vr, 
(AP-)+(z) = (M-) + (z) M, Z £ (-oo,0), 



(1.19) 



with the same asymptotics as z^oo. 

Note that the problems coincide at a = 0. In either problem, (M i ~)±(z) indicate the limits of 
M^{z) as z approaches either £"* or from the positive or negative sides (the precise sense in 
which the limit holds is discussed later). Last, 03 denotes the third Pauli matrix, (J _ ( j ) ). 

A large part of this paper is dedicated to the proof that there exist unique solutions to RHP'* 
and RHP^ . Granting that, we may now define the functions comprising the limiting kernels (|1.15p 
and (|1.16p . hereafter denoted A a (x,y) and M a (x,y). 

Definition 1.2. For — §7r < argz < |-7r and z ^ [0, a], 

(/;W>?:W) = iefe-^ 3/2 ((M-)n(z), (M-) 21 (*)) . (1.20) 

V 27T 



Similarly, 

(/£(*),£(*)) = -^e^e-^ 3/2 ^ 1/2 )(( M -)n(z),(M-)2i(z)) (1.21) 



/2tt 

/or a// z — |7r < argz < |-7r. 
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I) 



Figure 1: The contours and S^" for RHP"' and RHP^ 

One concludes that f^( x ) an d ga( x ) are rea l analytic for x > a and x > 0; the diagonals 
A a (x,x) and B a (a;,x) for x > a are therefore well defined. As for their behavior as functions of a: 

Theorem 1.3. Both A a (x,y) andM a (x,y) are continuous functions of a for fixed x, y. Continuity 
holds down (or up) to a = from either side. 

Finally, while we have not expressed the limit kernel in terms of known special functions, we 
do have the following asymptotics. 

Theorem 1.4. Uniformly for z in compact sets of (0, oo), 

f:(z) = Ai(*)(l + 0(e-* 3/2 )), £(z) = Ai(*)(l + 0(e~ a3/2 )) (1.22) 
as a— > + oo ; while 

£■(*) = (lal-I^JoC^dal-^^ + Odar 1 )), (1.23) 

i(^) = ^(^(H-^^+odar 1 )) 

as a— > — oo. Io(-) is the modified Bessel function of the first kind. 

After describing applications of Theorem 11.11 to the limiting distributions of the largest eigen- 
values for Unitary ensembles and some possible extensions, the analysis begins in Section 2 where 
the RHP connected to the polynomials {pk, n } is introduced. Section 3 subjects this RHP to a 
series of transformations, following the Deift-Zhou method of steepest descent [6j. A local analysis 
for the problem in the vicinity of z = 1 in terms of RHP^ and RHP^ is detailed in Section 4. 
With these parametrices, Theorem II .11 is proved in Section 5. Section 6 is devoted to the existence 
question for RHP^ and RHP~*; this is accomplished by a general vanishing lemma argument. Sec- 
tion 7 establishes several properties of those solutions, including their continuity and asymptotics 
(Theorems 11.31 and ll.4p . 

Remark The results here should be compared with those in the recent paper [3] which considers 
the following set-up: Take V regular with right-most edge of fj,y placed at the origin, and seek the 
asymptotics of the corresponding orthogonal polynomial kernel for the weight e - ™^^) restricted to 
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(—oo,0]. This is the same starting point as our problem. In [4j though, a parameter c = c(n) > 
is introduced in the weight as in cV which, when adjusted, can move the edge of the support of 
\iy away from the origin or push it in, creating a "hard-edge", or square-root singularity at the 
origin. While different at all finite n, this device is qualitatively the same as our choice of the sign 
of a, and should lead to the same phase transition in the limit. By a quadratic transformation, 
the authors of [1] are able to write the limiting kernel for points x and y to the left of the origin 
in terms of Painleve II. For the probabilistic motivations here, it is the kernel to the right of the 
critical point which is important. By analogy this should correspond to the kernel in [4] along the 
imaginary axis; the precise relationship remains to be worked out. 



1.1 Janossy densities and the distribution of the largest eigenvalues 

Denote the ordered eigenvalues of M by Ai > A2 > • • • . The well known gap formula for determi- 
nantal ensembles with kernel K n states that: for any BcK, 

P I there are exactly m eigenvalues in B ) = — — — det( I — 9K n ls I . (1-24) 

\ / ml da m \ ) e=i 

With m = 0, this formula together with the limiting result (|1.13p impies 

Urn p(Ai < 1 + — =det(l-^r« / l [a ,oo)), (1-25) 

with KAiry standing in for the L 2 -operator with Airy kernel, see again [10] for a full proof in the 
case of polynomial V. The celebrated result of Tracy and Widom ([20], with extensions in |21j ) 
provides a closed form for this Fredholm determinant, to wit, 

det(l - K Airy l [a oo ^j = exp(- J (s - a)u 2 (s)ds^j = F TW (a), (1.26) 

in which u(s) is the unique solution of Painleve II with u(s) ~ Ai(s) as s— > + 00. 

Formulas for the limiting distributions of the scaled A2, A3, etc, also exist. This is also found 
in [2D], though note there the asymptotics are only taken on for GUE. To explain, first replace 
the appearance of u(s) in (jl.26p with u(s;6) determined by the same equation but with u(s;6) ~ 
\/#Ai(s) at infinity, and denote the corresponding exponential function F(a;9). Then, following 
(|1.24p . x da F(a; 6) evaluated at 9 = 1 yields the (limiting) probability of there being exactly 
m eigenvalues larger than a. The corresponding distribution functions can then be built in the 
obvious manner. 

The Janossy densities provide a different path to the law of the scaled largest eigenvalues. From 
the definition (jl.5p . we have that 



dec 1 dx ■ 

Kt,k<m 



P ^exactly m eigenvalues in B^j = P^no eigenvalues in B 

X ~T / ••" / det L n,B(x£,X k ) 

m\ J B J b L 

and, assuming (jl.25p . one can obtain the following from Theorem II. li 
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Corollary 1.5. With M a equal to A a for a > and B a for a < 0, 

m— 1 



lim P(A m < 1 + ■ 



cyn 



2/3 



F TW (a)x^2— ■■■ det[M a (x e ,x k ))dx 1 ---dx n . (1.27) 
n =o n ' Ja Ja 

This describes the general limit distribution as that of the largest eigenvalue modulated by a 
finite sum of (standard) determinants with a universal kernel. Again, at this point the kernel is 
only defined in terms of a pair of RHPs, and so the above form of limit law is far from optimal. 

Proof of Corollary \1.5[ We provide just a sketch, using estimates developed below. To produce the 
second factor in (|1.27|) . one must pass the point- wise convergence of the kernel established Theorem 
11.11 under the integral. The fast decay of the exponential weight will control the integral at infinity, 
while more information is needed to deal with the integral near a. In particular, here one wants 
to show that n~ 2 / 3 L n ^ a (l + n~ 2 / 3 x, 1 + n~ 2 l 3 x) is uniformly integrable over x S [a, a + e]. Local 
forms of the solution of RHP^, provided in (|6,8p for a > and (|6.1ip - (|6.13p for a < 0, show 
that the first columns of along with their derivatives are bounded down to a or the origin 
(from the right). By Definition 11.21 we then see that, even on diagonal, the limit kernel (M a (x,x) 
= f' a ( x )9a( x ) ~ fa(x)g' a (x), neglecting the (<-) -superscripts) is integrable near a. Next, (|5.3p and 
(|5.5p express the finite n kernel in terms of certain (well-behaved) auxiliary functions and the first 
column of M^(z). A simple analysis of those auxiliary functions shows that n~ 2 ^ 3 L nyCt inherits the 
integrability of M*~ and yields the result. □ 

Improved asymptotics for the Janossy kernel might also provide estimates on the speed of 
convergence to the Tracy- Widom law (analogues of either the Berry-Esseen estimates or Edgeworth 
expansions for the classical central limit theorem). Results of this type are important in multivariate 
statistics, and have already been established for GUE and the related LUE in [T3] and [3]. The 
case of unitary ensembles with non-quadratic potentials has not been explored. Note however from 
the general formula (jl.6p we have 

d . _ /. a \ d 



— lo g ^^A 1 < 1 + — = — logdet^-A^l^^^^j (1.28) 

1 / a a 

7-L<n,a I 1 H TTfT: > 1 ~\~ 



cyn 2 /^ ' V cyn 2 / 3 ' cyn 2 / 3 ' 

A similar expression at n = oo is a first step in the derivation of (|1.26|) . and the limiting kernels 
(|1.15p and (|1.16p are not surprisingly tied to the resolvent kernel of the Airy operator, see [9j and 
[22]. More to the point, a suitable expansion in n in (jl.28p would bound the convergence speed. 
From Theorem 1 1 . 1 1 one anticipates the rate is n~ 2 / 3 , and that is just what is proven for GUE and 
LUE. Of course, carrying out the suggested program requires sharp asymptotics of L n>a (x, y) along 
the diagonal {x = y = a). An estimate of the form 

1 _ / a a 



L n!a (l + ^—,l + ^—)-M a (a,a) <n~ 2 / 3 (P(a) 



cyn 2 / 3 \ cyn"i' J cyn 
with 4>[a) integrable at positive infinity would, for example, more than suffice. 
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2 First RHP and introduction to the calculation 



The starting point is the RHP characterization of orthogonal polynomials due to Fokas, Its and 
Kitaev |13j . Fix a half-line V = (— oo,c] and consider the polynomials 

{Pk,n = 7k,nX k + • • • } orthonormal with respect to w n {x) = e~ nV ^ for x G V. (2-1) 

Then, the RHP reads as follows. 

RHP for Y: Seek a 2 x 2 matrix valued function Y(z) = Y n (z) such that 

Y{z) analytic in C\T, 

Y + {z) = YAz){ 1 MZ) ), z€T, 

K J V ' \ 1 / (2.2) 



r(*) = (j + o(i)) ( 



z n o 
z" n 



Z— KX3. 



The second, or "jump", condition, is read as 

Y±(z) = lim Y((), z £T, z G C±, the upper or lower half-plane. (2-3) 

2'— >2 

This can be understood in the sense of continuous boundary values for z6T away from the endpoint 
z = c with the additional condition that 

/ 0(1) 0(log|z-c|) \ 
r (zj = , 2 near c. 

V o(io g |*-ci) y 

With that said, the basic result is that the unique solution of this RHP is given by 



Y(z) = | 7n ' n 7n ' n V J > x | , (2.4) 

-27ri 7 n _i in p n _i jn (z) -27ri 7 n _i, n C\p n -x, n w n j (z) 

where C denotes the Cauchy operator on F: 

Cf(z) = Cxf(z) = J-flW-ds, zi S, 
2m J-£ s — z 

for any contour EcC and function f(z) G £ 2 (£, |cfe|). 

Note that (|2.4p contains the (iu n , T) orthogonal polynomials of degrees n—1 and n in its first col- 
umn. It follows that the kernel of interest, L n ^ a , may be expressed entirely in terms (Yu(z), Y"2i( z )) 
where T is now a function of both n and a: we have in particular, 

r = T n , a = ( — oo, 1 + ac7}vT 2 l z 



(2.5) 

and will use the additional shorthand c ni0l = acy l n~ 2 ^ for the (moving) endpoint. 

The analysis of Y for n— >oo entails a series of transformations, Y i— ► T i— ► S 1 i— ► i?, in order 
to obtain a RHP for which is normalized at infinity (i.e., R(z)^>I as z^oo), and has jump 
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matrices which are uniformly close to the identity as n^oo. Afterwards, unfolding this series of 
transformations will produce the asymptotics of Y. We are primarily concerned with the behavior 
of Y(z) in the vicinity of z = 1, for which we will build local parametrices. The basic program is 
identical to that in the analysis of the RHP connected to orthogonal polynomials over the full line 
in [IT] or [12]. Novel here is that the problem will follow two different paths, depending on the sign 
of a. 

2.1 Equilibrium measures and the (/-function 

The first transformation, Y i— > T, rests on properties of the the density which minimizes the 
analogue of Iy in (|1.9|) . We begin by recalling several properties of ifty, the "unconstrained" 
equilibrium density connected to the analysis of the full-line orthogonal polynomials. 

As indicated in the introduction, the support of tpy is a union of (N + 1) disjoint intervals and 
we normalize the right endpoint to sit at 1. The intervals of support are referred to as the bands; 
the complementary N intervals making up the gaps. Following [11] , the interior of the support is 
denoted by 

N+l 

J = (J (h~i,ak), 
k=i 

and the density tpv can be written, 

Mz) = 7^-R+\ z ) h v{z), for z e J, (2.6) 

in which 

N+l 

R(z)= H(z-b k ^)(z-a k ), (2.7) 

k=l 

and hy is real analytic on M. Here, the branch of R(z) is chosen so that R(z) behaves like z N+1 as 
z^oo. 

With this, the first transformation of the RHP for the unconstrained polynomials is based on 
the introduction of the g-function, 

g(z) = / log (2; — x)ipy(x) dx. (2-8) 
Jr 

The function g{z) is analytic on C\(— 00, 1] and has the following properties. First, there a constant 
t so that 

g+(z)+g-(z)-V(z)-£ = 0, for z £ J, (2.9) 
g+{z) +g-{z) -V{z) -t < 0, for z £ M\J. (2.10) 

(The strict inequality in (|2.10p is our final regularity condition.) Second, it holds that 

g+(z) — g~(z) = 2ixi j ipy(x)dx, for z & (— 00, 1). (2-11) 
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That is, g+{z) — g~(z) is purely imaginary on R and constant in each of the gaps, with the more 
detailed picture being: 



2iri, z G ( oo, bo), 

2ni J hh i/i v (x)dx = 2mVL k , z G (a k , b k ), (2.12) 

0, z G (l,+oo). 



9+{z) ~ 9-(z) = < 



Appraisals (|2.9[) through (|2.12p are basic consequences of the Euler-Lagrange equations for (|1.9p . 
as is explained in Section 3.2 of [TT] . 

With {pfc,n} we are working on the n-dependent interval T n ^ a , but may proceed in a like manner. 
For each integer n and real a, the old reasoning will show that infimum of 

I n ,M= I lo S i Mx)dn{y) + / F(x)d/i(x) (2.13) 

is uniquely achieved. The minimizing density is however qualitatively different a > or a < 0, and 
is denoted by ijjy or tpy respectively. 

2.1.1 The case a > 

When a > 0, Tn,a contains the support of the full-line minimizer ipVi and we state without proof 
the following. 

Lemma 2.1. It holds that ipyi 2 ) = Vv^)- Thus, by assumption U.12\) . 

i; v (x) = (l-x) 1 / 2 p v (x), x€ {1-6,1), (2.14) 

for all small 5 > 0. Here, /3y(z) is analytic in a neighborhood of z = 1, 0v(l) > 0, and for later 
IPC remark c v = (/?y (l)/2) 2/3 . 

Note that a definition of Pv(z) is implicit in (|2.6|) . Also, (|2.14|) has been set apart as this regime 
is of central importance in what follows. 

2.1.2 The case a < 

For a < 0, one attains different minima in (|2.13p and (|1.9p . Still, using the assumed regularity we 
will show that, for all n large enough, the support of ipy consists of N + 1 intervals, 

N+l 

J n ,a = |J (b k -x(n,a),a k (n,ai)\, 

k=l 

with oat+i fixed at c njQ , = 1 + a/(cyn 2 / 3 ). Now setting, 

= {z ~ bN \ T\(z- b k ^)(z - a k ), (2.15) 
[z - a N+l ) AA 

(from here on we suppress the (n, a)-dependence of the endpoints) we have: 
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Lemma 2.2. For a < and V satisfying $1.12]) it holds, 

i ~ 1/2 r i 

VV( Z ) = 7T-( R n,a( z ))+ ( C ™,« ~ z)h Vn a (z) + C n Q , (2.16) 
2,711 L J 

on its support, with 

C n , a = \{^ m )h v {l) + 0{n-^), 

where hy is as in \2. 6\) and hv <n ,a( z ) *s rea >l analytic and tends to hy as n— >oo. The analogue of 
firm ) reads 

^y(x) = (c n>a - *) 1/2 /3S,» + I(-^)(c^ - xy^P^Jx), x € [1 - 5,Cn >a }, (2.17) 

where Py'^ a (z) are analytic in a neighborhood of z = 1 and lim n ^oo Py'^ a (cn,a) = ^v(l)- 

With i n hand, we define a ^-function exactly as in (12, 8D . The basic relations (12, 9h through 
(|2. 10p remain valid, though with adjusted values for £ and the S7's. 

Example (GUE). It is instructive to first spell out the computation for GUE, or when V(z) = 2z 2 . 
The unrestricted minimizer has one band of support, [—1,1], and is given by the semi-circle law, 
iI>gue(x) = — x 2 . Now fixing the right edge at a = 1 — e for small e > 0, there is still one 
band [6, a] with b and the adjusted density V'gc/.e to be identified. If we put 

G(z) = 1 / 



ds, 



s — z 

differentiating the relation (|2.9I) produces the scalar RHP: 

4i 

G + {z) +G-{z) = —z, z€[b,a], and G+(z) - G-(z) = 0, z e R\ [6, a] . (2.18) 

Introduce = R £ ( z ) = fEl an d multiply both sides of (|2.18p through by the square-root of this 
object. The RHP is then transformed into a standard form, and one finds 



c W = vgw r (Wj (2 . 19) 
2m • / » (/^)) + s - 2 



subject to the single moment condition, 



which holds since >i as z^oo. The integral (|2.20p is easily computed, and b = b(a) = 

|(a — 2V a 2 + 3) for positive a < 1. Also, by properties of the Stieltjes transform, it holds that 

^gue(x) = Re (G+(x)) = ^y^T^— 2~ ~ x ) ' x e (°' 6 )- 
Now, for e small and x near a, 

= ((1 - e - x) 1 ' 2 + |(1 - e - x)-V2) (2^2/vr) (1 + C( e )) f (2.21) 

after substituting a = 1 — e, b = — 1 + C(e 2 ) in the previous display. This provides a model for the 
general formula (12. 17ft . □ 
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Proof of Lemma \2.^A As in the previous example, set e = ac v 1 n 2 / 3 with oat+i = a7v+i(e) = 1 — £. 
Following the standard approach, the density is given by 



Vv, e (z)=Re((G £ ) + (z)). 



where, 



and 



<M0 - ± VW in I - «)— , (2-22) 

2vrz (z - ajv+i(e)) Jj e (y/R £ (s)) + s-z 



R e{z) = Yl ( z ~ ~ Ofe(e)), J e = IJ afc(e)). 

^=i fc=i 

We have normalized in this manner as, when e— >0, the endpoints Ofe(e)) converge to their 

positions in the full line equilibrium density t/Vj and R £ {z) and J e converge to R(z) and J defined 

in (|2.7j) and directly above. 

The integral in (|2.22p over J £ can be replaced by 

for C a clockwise oriented contour surrounding both J £ and z. Now, 

iV'(s)/2-K ds f iV'(s)/2n ds f iV'(s)/2n 
(ajv+i(e) - s) = (a N+ i(e) - z) / — / — 



= (a N+1 {e)-z)h v , £ (z) + C £ , (2.23) 

where hy >e (z) and C £ are the same as hv, n ,a( z ) an d C nja in the statement of the Lemma. Since 
hy(z) figuring in the definition of ipy is exactly lim £ ^o hy,e (recall (12, 6p ). it is left to prove that 

C £ =[^L^-U Z^L_Li! + o<f). (2.24) 

Indeed, 

by a moment condition for the full-line (or "free") problem, and the second integral in (|2.24p is 
exactly hy(l). Further, the endpoints {fefc(e)} and {afc(e)} turn out to be real analytic functions of 
e, and we have 

e f V'(s) a' k (0) , y fc -i(0) \ da 2 

The advertised 0(e) term will then arise from otv+i(0) = 1 and a' N+1 (0) = — 1, and additional fact 
that all other endpoints have vanishing first derivative at e = 0. 

The idea behind verifying these last claims is to view our constrained problem as a perturbation 
of the free problem at the modulation point P, where qtv+i = 1. Returning to (|2.22p set R(z) = 
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R(z)/(z — djv+i) 2 , suppressing the dependence on e > 0. For the constrained problem, the system 
of 2N + 1 modulation equations determining the endpoints are 

Tj = f ^=kl s j ds = 0, < j < N - 1, (2.26) 
f R(s) , 



T N = [ ^kLs N ds = ±, (2.27) 



the so-called moment conditions along with the integral conditions 



N k = G(z)dz = 0, 1 < k < N, (2.28) 

J 7k 

where 7& denotes the loop around the cut (&&_!, a&). Adding the additional relation 

T_, 5 / pfl d ' =0 (2.29) 

gives the system for the free problem. In particular, the moment conditions for the full problem 
read fj V'(s) / \J R{s)^ds = for < j < N and (tti)^ 1 for j = N + 1, and T_i takes you from 
the system (|2T2lJj) - ([2T2Tj) to this one. 

Our regularity assumption for the free problem implies that the absolute value of the (2iV + 2) x 
(27V + 2) Jacobian of the map {cij , b^} i— > {Tj , N^} is bounded below by a positive constant at point 
P. A full proof may be found in Section 12 of [IS]. In order to conclude all other endpoints are 
real-analytic functions of oat+i, and so e, we must show that the (2iV + 1) x (2N + 1) Jacobian of 
the modulation equations for the constrained problem is also bounded below by a positive constant 
at point P. (And then invoke the implicit function theorem). Clearly, this will hold if 

9Tj „ , „ . . , dN k 



for < j < N, and — = for 1 < k < N, 



da N+ i " da N+1 

as, if so, the (2A^+ 1) x (2N + 1) minor must be non-degenerate, less the full Jacobian is. Note also 
that by a simple application of the chain rule, this will imply the vanishing of the first derivatives 
of all endpoints other than aN+i at e = 0, a fact we used above. 

Finally, g a ^ +1 Tj = Jj V'{s)/ \J R{s) + s^ds and the latter vanishes for j < N by the moment 
conditions for the free problem. Also, by linearity the vanishing of the derivative of any will 
follow from (P) = 0. We compute 

dG _ -1 , f V'(s) ds | 1 f V'(s)(s-a N+1 ) ds 



da N+L Airi* Jjy/R{s) s-z Am (z - a N+l ) J j ^/r[ 



1 yR{z) f V'(s) 

-ds, 



Am (z - a N+ i) Jj y/R(s) 

and note the last integral vanishes at point P as the integral appearing reduces to that in condition 
(|2T29|) . □ 
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3 Steepest descent 



3.1 First transformation Y i— ► T 

As in Section 3.3, we define T{z) for either a > or a < by conjugation, 

T(z) = e~^^ Y(z) eir^e-nsW^ ( 3-1 ) 

1 



where g(z) is the log-transform of either tpy or ipy, and we recall that 03 = y ^ j. The jump 

matrix for Y{z) is transformed into, 

/ e -n(g+(z)-g-(z)) e n(g+(z)+g-(z)-V(z)-e) \ 

V ^=[ o c -*(M.)-*-M) J' Zer ^ 

Next, using the relations (|2.9p and (|2.12|) satisfied by g(z) and the fact that e n9 ^ ~ z n for z^oo, 
we find that T(z) is the unique solution of the following RHP. 

RHP for T: We seek analytic in C\r n a , with jump relations, 

/ e -n{g+{z)-g-(z)) j_ . 

T+(z) = T_(*)( ... ), zeJ 

T + (z)=T_(*)( q e2Mj ), ^€(0^6^ = 1,...,^ (3.2) 

j e -n(g+{z)+g-{z)-V{z)-l) 

T + (z) = T-(z){ ^ J, z < b or oat+1 < z < c njQ ,, 

and asymptotics, 

T(z) = I + o(^j z^oo. (3.3) 
That is, we now have an RHP normalized at oo. 

3.2 Second tranformation T i— > S 

This step is the descent, transforming the oscillatory diagonal entries of the jump matrices in 
the RHP for T(z) into exponentially decaying off-diagonal entries in an equivalent problem for a 
function S(z). 

For the a > case we follow [11] without change. For z £ C\r njCf in the region of analyticity 
of hv(z), recall (|2.6p . define 

0(z) = f R l/2 (s)h v (s)ds, 

J a N+ i 

where ajy+i = 1 and the path of integration does not cross T n ^ a . From (|2,9p and (|2.12p we have 
that, for each z £ (bj—i,a,j) C J, 

ra-N+l 

g + {z)-g_{z) = 2ni ^ v {s)ds (3.4) 



R 1 * 2 (s)h(s)ds + 2m / V(s)ds 
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That is, —4>(z) and 4>(z) are analytic continuations of g+{z) — g~{z) above and below each band. 
Also, 0+ and 0_ are purely imaginary on each band and an easy exercise using the Cauchy-Riemann 
conditions shows that, 

Re <p{z) < 0, for small Im (z) ± and Re (z) £ J. (3.5) 

Since the ^-function tied to ipy (for a < 0) satisfies the same basic relations we can define extensions 
for g+(z) — g~{z) = 2iri f^ n ' a t^y(s)ds, and so (j), in the same way. 

With these properties of eft in mind, the jump contour is deformed off the line by opening a lens 
around each band based on the factorization 

/ e -n(9+(z)- g -(z)) 1 \ / 1 q \ / o 1 \ / 1 q\ 

I o e n (g+{ z )-g-( z )) J I e n<t>-(z) i J \ _i oil e n< t'+^ 1 / 

= BZ l {z) ([ 3 i X \ B + (z). 

Set, 

I T(z), z in the exterior of each lens, 

S(z) = < , (3.6) 

I T(z)B± (z), z in the upper/lower part of each lens, 

in which what is meant by a lens, and the resulting contour with lenses S nQ! , is spelled out in 
Figure [2J We are led to: 

RHP for S: S(z) is analytic in C\£ nja , satisfies S(z) = I + 0(~) as z^oo with z £ S raja , along 
with the jump relations, 

S+(z)=S-(z)\ ), zGS nia nC±, 



S+Oz) = I ] , z G J, 

e -27ri^ e -n(fl+(*)+ fl _(*)-V(*)-fl . ( 3,7 ' ) 

q e27r . n . J, ze(a j ,b j ),j = l,...,N, 

, 1 e -n(s + (z)+ fl -(z)-^(«)-^) . 
S+(z)=S-(z)( | J, z < 6 or aAT+i < z < c njCf . 

Note that (|3.5p implies the factor e"^W in (|3,7p decays exponentially as n^oo. Further, in the 
regular case V the exponent g+ (z) + g- (z) — V(z) — £ appearing in the third and fourth jump matrix 
is strictly negative in the gaps or past the ends of support. It follows that the corresponding entries 
also decay exponentially as n^oo. 

3.3 Model Problem 

From the discussion at the end of the previous section we expect the leading order asymptotics to 
be governed by the 2x2 matrix P°°(z) which solves the following model problem. 
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Figure 2: The new contour E njQ; (pictured for a > 0) with a lens opened around each band. 

RHP for P°°: P°°(z) is analytic in C\[b ,a N+1 ], P°°(z) = I + 0{l/z), as z^oo, and 

, , ( 1 \ 




o \ (3 8> 

Pf!:)-P-!--)| Q e _ 2 , jnj J, *£(«,,&,), J = 1,...,JV. 

Though we are led to this problem by considering the n— >oo for the jumps of S'(z), the bands and 
gaps over which the jumps of P°° are defined should still be taken in their finite n positions for the 
a < case. 

While the particulars of P°° will not affect the parametrix we eventually construct about z = 1, 
it is required to demonstrate that the above problem does indeed have a solution. Fortunately, this 
has already been accomplished in [TT] . where it is proved that (|3.8p has a unique solution satisfying 
detP°°(z) = 1. 

3.4 Last transformation S i— > i? 

Since the convergence of the jumps for £(£) to those for P°° is not uniform near the endpoints, we 
have to perform a local analysis at each of the endpoints a,j, bj. The analysis at oat+i is particular 
to the present endeavor and is the subject of the next section, for the rest though we may again 
refer to [TT] . 

For n large enough, each xq = a,j,bj is regular and there will be no interior "singular" points. 
Surround each xq by a small disk and consider the set of local RHP's: 

P xo (z) analytic in {\z — xq\ < e'}\T, ntOC for a e' > e, 

P xo (z) and S(z) share jump conditions on T> n ,a H {\z — xo\ < e}, (3-9) 
P a;0 (z)(P oo )~ 1 (z) = I + 0(n~ K ), uniformly for \z - xq\ = £ with a K > 0. 

The last condition matches asymptotics of P Xo (z) to those outside its disk. Granted a solution we 
define 

i S(z)(P 00 )- 1 (z), z outside the disks, 

R(z) = < ):,_,, ' (3-10) 

I 5(z)(P.) (z) z inside the disks, 

in which P,{z) stands in for whichever P Xo (z) corresponds to the given disk. At all xq ^ ajv+i; it 
is well known P»(z) is given explicitly in terms of Airy functions; the connected RHP is in fact 
RHP~* with a = oo, and the error exponent is re = 1. Assuming that parametrices for z = ajv+i 
exist, R(z) will be analytic off the system of contours described by Figure [3l (While there is an 
isolated singularity in the second column of R{z) at a^+i = c n ,u traced back to that in Y(z) at 
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the same point, it is logarithmic and so removable). Next, since S(z) and P°°(z) are normalized 
at infinity and det P°°(z) = 1, it follows from its definition that R(z) is also normalized at infinity. 
These facts are key ingredients of the result that, uniformly for z £ C\£#, 

R(z) = 1 + 0(n- K ), n^oo, (3.11) 

with proof identical to that in [IT]. It follows that R(z), and also -^R(z), are uniformly bounded 
for large n, and from (|3.10p that det(R(z)) = 1. 



Figure 3: Part of the contour 



4 Paramatrices at the right edge 

At last we build paramatrices for the RHP's (a > or a < 0) in a neighborhood of z = a^+i as 
n— >oo; these are described in terms of the solutions of RHP^ and RHP^ . 

When a > 0, we have that cln+i = 1 < c n , a . It is convenient at this point to bring in the fact 
that, 

g+(z) + <?_(» + V(z) - £ = -J" R 1 ' 2 {s)h v {s)ds = -<j>(z), (4.1) 

for real z with \z\ > 1. Then, for a fixed U £ = {z : \z — 1| < e} we require a 2 x 2 P~^(z) which is 
analytic in U £ \T< nta and satisfies, 

(J~(*)) 

with 

P^(z)(P°°(2))- 1 =I + C»(i), ^^\S n , a . (4.3) 

The asymptotics (|4.3|) entail a matching condition between the inner and outer solutions. 

If instead a < 0, the endpoint depends on n as in a/vr+i = 1 + a/cyn 2 / 3 = c n;Q . It is convenient 
though to keep the same neighborhood U £ fixed about z = 1 and the problem is to find P*~(z), 



(P-(*))- 
(^(*))- 



1 > 

e ncf>(z) x ^ 

i\ 

-i y 

1 e -n(t>(z) 
1 



z g (E„, a n u e ) n c±, 

z e (1 -e,l), 
ze [1,1 + 



(4.2) 



273 7) 
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analytic in U £ \E na with jump conditions 

(P-(z))+ = 

and again, 



oo / ^\ \ — 1 



J"-(z)(P°°(*)) 



i + o(- 

,n 



1 

e n<j){z) j 

1 \ 

-i ; 
i 



(4.4) 



(4.5) 



These problems are now mapped onto those for Af^(£) and M*~(Q. For z G U £ define two 
changes of variables z^Q = Cn O 2 )) via 



and 



Also, set 



n<j)(z) = — C 3//2 , for a > 0, 



n 0(z) = -C 3/2 + 2aC 1/2 , for a < 0. 



^(z) = P°°(z)- 



1 



(C»(*)) 



a-3/4 



2i V 1 1 

Lemma 4.1. T/ie i?#Ps fggj, g^PJ and are solved by 

P-(z) = E n (z)M-(Q(z)) and P^{z) = E n (z)M*~(Q(z)). 



(4.6) 
(4.7) 
(4.8) 

(4.9) 



Here Cn( z ) and Cn( z ) is given by \4.€ty or \4. 7| ) respectively. Also, M^{Q) represents the solution 
of RHP^ in which the jump over [0, a) is replaced by the same jump over [0, a n ) with a n ^ a as 
n— >oo defined below in ft4.1!fy . 



Proof. For a > 0, from (|2.14p we have that: with z close to 1, 

n(j}(z) = n [ {s- l) 1/2 (3 v {s)ds = \n{s - 1) 3/2 (3{z), 

where (3{z) inherits the analticity properties of (3y. Now set 

( n (z)=n 2 / 3 (z-l)(P(z)) 2/3 . 



(4.10) 



(4.11) 



Along with being analytic in U £ , /3(1) > 0, and the branch may be chosen so that, (/3(1)) 2/3 > 0. 
Then, by choice of e, C, n {z) maps U £ one-to-one and onto an open neighborhood of £ = Cn(l) = 0. 
Further, ( n (U e flR) CM, Cn(U £ n C±) C C±, and those parts of the z-contour (E n>a n U £ ) n C± can 
be chosen so that their images are arg£ = ±§7r. Last, z € [1, 1 + a/cyn 2 / 3 ] is mapped to 



0,a(/?(c n , Q )) 2/3 /cy)l = [0,a n }^[0,a((3v(l)/2) 2/3 /cv] = [0, a] 



(4.12) 



as n— >oo, providing the definition of a n . 
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When a < 0, we have from (|2.17p that, 

™Kz) = nf ((^-g^t^^^-g-^wjd. (4.13) 
= {^n(z - c n , a f/* + + Q{n-^)\ (z - c n , a )^}p{z). 

Note the change of sign: (s — Cn,a)+ = — (s — Cn, a )+ ■ Here again (3{z) is analytic and positive 
in a fixed neighborhood of 1, and so c na , for all large enoug h n. (That f3 { y^ a (z) are analytic near 
1 and differ by 0(n~ 2 / 3 ) at z — c n Q , is used.) Choosing the 2/3-root of (3 positive in that same 
neighborhood, it follows that (|4.7p has the required properties: U e is mapped one-to-one and onto 
a neighborhood of C = Cn(c n ,a) = 0, C«(^ H f/ e ) is real, and the segments of £ n)0 can be chosen to 
map onto arg£ = ±§7r. 

Plainly, M^{Q n {z)) and M*~(£ n (z)) satisfy the jumps specified for P~*(z) and P^(z). Next, as 
fully explained in [12], E n {z) is analytic in U E and so does not affect the jump relations. Briefly, 
P°°(z) has a singularity of the form (z — aN+i)~ az ^ 4: at the right edge, and this is compensated 
by the appearance of {Cn{z)Y 3 ^ in (|4.8|) and the fact that (, n (z) Cra 2 / 3 (z — ajv+i) whether a is 
positive or negative. The asymptotics (|4.3p and ()4.5p follow from the fact that there is a constant 
c > with |Cn(^)| > cn 2 / 3 uniformly on \z\ = e and n— >oo and the behavior as rj— >oo of M*~(Q 
stated in (fTTTHD . □ 

Remark 4.2. Along with E n {z) being analytic in a neighborhood of z = 1, it also follows from the 
form of P°°(z) that both E n (z) and 4-E n (z) are uniformly bounded in that neighborhood. Further, 
det E n (z) = 1 since det P°°(z) = 1. 



5 Proof of Theorem 11.11 

The derivation below borrows heavily from [TTJ. First note the expression for L n ^ a in terms of the 
Y matrix (defined in (|2.2p ): 

L na (z,w) = - J- e -^W-|v(«) ^ii(^2iH-l2i(^nH^ 
2iri z — w 

We are interested in the asymptotics of the above kernel for real z and w to the right of ajv+i, the 
endpoint of support of the density of states. Unravelling the sequence of transformations leading 
from Y to R, we have by (ED, jMJ and (l3TT0|) . 

Y(z) = Y na (z) = V2^e^e^ 3 R(z)E n (z)M n>a (z)e^ z > 3 e-f a3 e n9 ^' T3 (5.2) 



in which we have made the definition: 

( pTi/i 



M^(Cn( z )) e '^ {z)cT3 for a > 



M nQ (z) = ^ i 2 f 4 . (5.3) 

n,a V 7 | ^ M -( C -( z )) e -|n0( 2 ) CT3 f or a < 
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Here Cn( z ) anc ^ (n are given by (|4.6p and (|4.7p , and is as in Lemma liTTl It follows that§] 

(5.4) 



Yu(z) 



I (Mn, a )2l(Cn(^)) 



Now, for either the (->) or (<-) case and real z > oat+i = 1 or oat+i = 1 — a/cyn 2 ' 3 , we have that, 

" \v(z) - l -= -(<?+(*) + g.(z) - V(z) -£) = -~<f>(z), 
recall (|4.1|) . Further, set 

K{z) = R(z)E n {z), 

and recall that both R{z) and E n {z) have determinant = 1, are analytic in a neighborhood of 
z = 1, and, along with their derivatives, are uniformly bounded there. Obviously, K{z) inherits 
these properties. Next change variables as in 

x . y 



x r . 



1 + 



cyn 



2/3 



, w y n = 1 + 



2/3' 



with x, y > a, and summarizing the steps thus far we have: 



1 



Cyn 



2/3 



Ln,a{Xm Vn) 



(5.5) 



1 



2m{x — y) 
1 



det 



det 



K(x r . 



e f e -§nV(xn)y 21 ( Xn ) e 2 re-^ y ^™)y 2 i(y ri ) 



(Mn, a )llW 

(M n , Q ) 2 i(x„) 



(M„, )ii(y„ 
(M_ a ) 2 i(jfo 



(x - y) 

Before further manipulations we record the following two facts. 
Claim 5.1. As n— >oo, 

C(x n ) = x(l + 0(n- 2 / 3 )), and Q(x n ) = (x - a)(l + 0(n~ 2 / 3 )), 

uniformly for x in compact sets of (a, oo). Also, 

f |^(l + 0( n -2/3 ))) a>Q) 
n<^nj - ncp nA x n ) - | ( ^ _ a)3/2 + ^ _ a)1/2)(i + 0(n _ 2/3)) q < Q) 

with the same uniformity in x as n— >oo. 
Claim 5.2. For all real x > a, 

MZ(x) - M~*(x) = 0(n- 2/3 ), 

and 

^(M-(x)-M-(x)) = 0(n- 2 / 3 ), 
as n^oo. The estimates are uniform for x a positive distance from a. 

1 



(5.6) 



(5.7) 



(5.8) 
(5.9) 



If z, w are taken to the left of ajv+i, there is an additional factor of 







-ncj>(z) 



in (|5.2p arising from opening 



the lenses in the T i— > S step. In that case, Yn and Y21 are linear combinations of the first and second rows of M„ i0 
respectively. 
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Proofs. The estimate (|5.7p follows directly from substituting the definition of x n into (|4.10p and 
(g33D : 

iuj>{x n ) = ^x 3 ' 2 x p(x n )/(c v f/ 2 , 



n<j>{x n ) = [-{x - af 2 + (2a + 0( n - 4 / 3 ))(x - a) 1 ' 2 } x ${x n ) / (c v ) 3 ' 2 , 



for a > 0, and, 

4 

.3' 

for ct < 0. Since f3(z) and /?(^) are analytic and = Cy > at z = 1, each of the rightmost factor 
above may be expanded as in 1 + 0(n~ 2 / 3 ) + ■■■. The same considerations lead to (|5.6p . 

As for (|5.8p and (|5.9p . Lemma 1731 below proves that M~*(z) = M^(z;a) is continuous in a for 
z supported away from £~*. The estimate (|7,4p obtained in its proof provides 

M~*(z; a) - M^(z; p) = O(o - /?), 

with the same holding for the derivative of the left hand side. Now, since M^(z) = M~*(z; a n ) and 
(|4.12p shows that a n = a + C(n~ 2//3 ), the claim is proven. □ 

Picking up the calculation, the matrix within the last determinant of (15. 5p is now written as 



K(x n ) 



(M n , a )n(a; n ) (M niQ ,)n(y n 
(M ni a) 2 i(x n ) (M n<a )2i(y n 



(5.10) 

i,a)2l[yn) J 

+K{y n )-\K{x n )-K{y n )) 



(M nia )u(x n ) 
(M n , a ) 2 i(x n ) 



Since the analytic matrix function K(z) satisfies det K(z) = 1, we see that the form of desired 
limit resides in the first term of (|5.10p . As for the second term, first note that since K(z) and its 
derivative are uniformly bounded for \z — 1| < e with a small enough e > 0, it follows that K(z)^ 1 
is bounded in kind and that K(x n ) — K(y n ) = 0(\x n — y n \) = 0(\x — y\n~ 2 / 3 ), by the mean-value 
theorem. Therefore, 



K(y n )-\K{x n ) - K(y n )) 



(M nia ) n (x n ) 0\f 0(\x - y\n~ 2 / 3 ) 
(M n , Q ) 2 i(x n ) i I 0(\x - y \n- 2 ' 3 ) 



Along with the estimates on the pre- factor, we are in the domain of analyticity of and M*~, 
which coupled with (|5.6p through (|5.8p . implies that M niQ (x n ) is uniformly bounded. 
The kernel now reads, 

j { , 1 / (M n , a ) u (x n ) + Q(\x - y\n~ 2 / 3 ) (M„,«)u(y n ) \ 

L n ,a{ X n,yn) = 7 r det (5.11) 

(x-y) \ (M n>a ) 21 (x n ) + 0(|x-y|n- 2 / 3 ) (M n ,«) 2 i(y„) / 
1 -detf (Mfl,a)ll(iCn) (Mn ' o)ll(2/n) V°(^" 2/3 )' 



(s-y) V (M niQ ) 2 i(x n ) (M niQ ) 2 i(y r , 
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for x, y bounded in (a, oo) and all large n. Instead of expanding the first term of the right hand 
side entry-wise, the second column of that matrix is subtracted by the first, 

L na ( Xn ,y n ) = T ±^det( ( M nMxn) ~ (M nM y n ) (M n , a ) u (y n ) \ 2/3 
' {x-v) \(M nta ) 21 (x n )- (M n>a ) 21 (y n ) (M n My n ) J 

which will allow an estimate uniform in x and y even as \x — y\ { 0. 

Using (|5.7|) and (|5.6|) and the analyticity of M < ~(z) (Rez > 0), one can check that: with 



O 



if! 2 / 3 ' 

with 6{x — a) = (4/3) (a: — a) 3 / 2 + 2a(x — a) 1//2 and • = 1, 2. For a > 0, the same basic reasoning 
gives: 

^((M n Mx n ) - (AC).i(*Ke* fl( * ) ) = O(^) 
with #(x) = (4/3)x 3//2 and again • = 1 or 2. It follows that 

(M^Uixn) - {M n , a ). x {y n ) = i ((M-).!(x - a)e¥^-«) - {M^). x {y - a)e^~ a ^ 

+0(\x- y\n~ 2 ^). (5.13) 

for a < 0, with the analogous statement for a > 0. Detailing how f)5. 13j) is employed in (|5.12p we 
recall the definition of (/£ ,g&) from (|1.2ip and write, 

Ln,a<oi.^ni Un) 

1 det / (/«r(x-a)-/-(y-a) + 0(|x-y|n- 2 / 3 ) fe(y - a) + 0(n" 2 / 3 ) ^ 



V 5a (a; -a) -9a(y~ a ) + - y|n~ 2 / 3 ) ^(y - a) + 0(n" 2 / 3 ) / 

TO / n, _J_ A J fa (x-a) -fa(y~a) 0(n~ 2 / 3 ) \ ^ 2/3 

= JB Q x,y ) + - -det , +0 ?! ' ). (5.14) 

\ X -V) \ 9a(x-a)-g^(y-a) 0{n 2 / 3 ) J 

The first term of the right hand side is the advertised limit kernel. To show that the second term 
is 0(n -2 / 3 ) uniformly in x,y in bounded sets of (a, oo) note that 

fa (g -a)- fg(y ~ a) _ Q ^ 
x-y 

for all such x and y since x t— » (M*~)n(x — a)e _ 2^ a:_Qi ) is smooth for x > a. The same is true if 
f£ is replaced by g^. This completes the proof for a < 0. 

In the case a > following the steps behind (|5.14p produces 

Next repeat the procedure: subtracting the second column from the first in the above determinant 
and now employing the estimates of Claim [5\2l will allow each appearance of (M^)-i to be replaced 
with the corresponding (M^).i with an overall 0(n~ 2 / 3 ) error. 
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6 Existence for the local problems RHP^ 

The general theory connects the the construction of a solution to a given RHP to that of a certain 
singular integral operator. In particular, consider the RHP (E,v): 

m(£) analytic in C\E, 

m+(C) = m_(0«(0. CSS, (6.1) 
m(C) = J + o(j), C^oo,C^S, 

in which v is continuous on E away from points of self-intersection and v{C,)—>I as £^oo along E. 
Next define the integral operator 

a/(o = c_ (/(«-/)) 

on L 2 (S, |c?C||), where C_, and C_|_ are the i-limits of the Cauchy operator: 

(C±/)(C)= lim ^- / -^ds, CGS. 
c'e± side 01 s 

With v G / + L 2 (E), C v is bounded from L 2 (E)^L 2 (E), and an I? solution to ([63]) can be 
construction out of a solution fi v £ I + L 2 (E) of 

(i - c^k = i (6-2) 

via 

m(C) = J + ^ / M*)M«)-/) da _ 
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For details behind these facts, |23] is recommended. Existence for the RHP (E, f ) would then 
follow from showing that I — C v is a bijection in L 2 (E). 

To apply this strategy to either RHP^ or RHP'" requires a preliminary step: it is not the 
case that M*~, or their corresponding jump matrices V*~, are normalized to the identity at infinity 
Therefore, we bring in 

-w-^M-i !)'"*•" (6 - 4> 

the fundamental solution of the twist problem: 

m+(C)=m_(C)( ^ M, forCGM_, 



and m(C) analytic in C\R_. With this, we define 



mo = \ M : {0m ' Ha for|cl>fi ' (m 

I M-(C), for |C| < fl, 
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with a fixed large and positive R. The contours for the pair of RHPs for appear in Figure [H 
Obviously, in each case the jump along the negative real axis has been removed far out. Further, 
along £2 and £4 and |£| > R, the new jumps 



V(C) = m(C) 




(m(z))-\ 



(6.6) 



with 9(C) = |C 2,/3 or 9(() = f £ 2 ^ 3 - 2 |a|£ 1//2 , still decay exponentially fast to the identity as (—>oo. 



Last, along the introduced contour |£| 



R, both problems have the uniformly bounded jump 




ICI=R 



\C\ = R 




Figure 4: The contours for (Y,^,V*~). 

We now have a pair of problems which fit into the above program (jumps are G L°° n (/ + L 2 )). 
The proof of existence now comes in three steps: to show I — is Fredholm, has zero index, and 
then that ker(I — C~^) = 0. This last point is established through a vanishing lemma similar in 
spirit to [IT] . Section 5. 

6.1 Fredholmness 

Fredholmness is implied by the following continuity condition holding throughout the contour. 
Moving clockwise about a point p on S, at which segments of the contour Si through with 
jumps v\ through meet, continuity at p is equivalent to 

I = v l {p) ±l v 2 {p) ±l ---v k {p) ± \ (6.7) 

in which the sign (±1) in the exponent is determined by whether the given contour points into, or 
out of, p. Additionally, this assessment is invariant of conjugations or deformations, see [23J. 

6.1.1 Criteria ACT for RHP^ 

The conjugation by m(C) in the exterior of a large disk not only produced a problem with decay, 
from the point of view of Fredholmness, it removed the discontinuity point at £ = 00. While new 
point of intersection £ = pi,P2 and P3 have been introduced to resulting contour by this move, each 
corresponded to a point of continuity for the (£*~, V*~) and thus remain so by the discussion above. 
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The worrisome points left on T," are then the origin and £ = a. At the origin, the jump matrices 
through S4 satisfy, 

-1 

/ 1 1 \ 

I 



1 1 
1 



1 
1 1 



1 
-1 



1 

1 1 



4 c 3/2 



and therefore (|6.7p holds. At the point £ = a, fix a small e > (e < a) and consider the local 
problem 

' 1 e _ 3 l 
1 

I (P a (0)- for CG (a,a + e] 

The jump being upper-triangular allows us to write down an explicit solution: 



for C G [a — e, a] 



4,3/2 

1 1 p 

2-7TJ JO- " 





-e s— C 
1 



(6. 



holding in L 2 and in the sense of continuous boundary values away from £ = a — e and £ = a. 
(Recall, C+ — C_ = I.) Next choose a positive e' < e and define 

M a (() = M^(()P a (() for |C -«!<£', 

leaving M~*(£) unchanged in the exterior of this disk. The effect of conjugating out the local 
solution is a new RHP with contour depicted below in Figure The point of discontinuity £ = a 
has been removed, with the introduced point of self-intersection at p = again automatically a 
continuity point, having arose from such by way of a conjugation. 




Figure 5: The introduced jump along the boundary of the disk D a ^ £ i is P a . 



6.1.2 Criteria (jSTfji for RHP" 

From the previous discussion for RHP" it is plain that we must only deal with the point Q = 
where the three ray of the contour T" come together. This is again handled by conjugating out a 
local solution local solution -Pq(C) of the RHP with jump V" restricted to Y>"C\{U £ = {Q : \(\ < e}). 
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Given P we set M o (0 = M^(Q P ((), for ( G C4' with e' < e, and M (C) = M<"(C) for C G C\i7 e /. 
The .R-ffP for M will satisfy ((6771), and it will follow that (e*-,V*~) is Fredholm. 
To construct Po((), consider 

M-(C)e-(f^ 3/2 +^ 1/2 ) CT3 = M^(C)e-^> 3 , C G C/ e , 

which has the constant jumps, 



o(o i r iii 



(6.9) 



Extending the jump contours on the right to infinity, we obtain the problem: find some Q((), which 
satisfies 



Q+(C) = Q-(C)| _° x Jj, C g 

CMO = ^-(U lli], argC = ±|7r, 



(6.10) 



-3' 

with Q(C) otherwise analytic. From [16], Section 6, we have: 

Proposition 6.1. With H^^-) and H^\-) denoting the Hankel functions of the first and second 
kind, and Iq(-) and Kq(-) the usual modified Bessel functions, the (un-normalized) RHP W.19\l has 
the solution: 

Q(C)=( l0( ^ V (6.11) 



for — §7r < argC < §7r, 



/or Svr < arg£ < 7r, and 



-vrVC(^ iJ )(v^C) vrVC(^)(v^C) 
/or — 7r < argC < — |vr. 

It follows that we can set -Po(C) = Q(C) e 5 ^^" 3 to form the needed local solution in \(\ < e and 
complete the proof. 



6.2 Index zero 

One consequence of the Gohberg-Krein theory of factorization of matrix- valued functions, is that, 
given Fredholmness, the index of I — Cy equals the winding number of det V over the contour, see 
[18] . But with V equal to either or V"*, det V = 1 and so that winding number is zero. 
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6.3 Vanishing Lemma 

Finally we show that ker(I — C v ^) = 0, first reverting back to the problem(s) tamed at infinity 
(IT, V* 1 ), recall (jfiSJ) and fl&i]). Suppose that ker(I - C~-) ^ 0. Then there exists ffi G L 2 (S) 
which satisfy 

(I-C^)/4 = 0, 

and so 

JTr s — C 

are L 2 -solutions of the RHPs: 

M,f (C) analyltic in C\XT or C\£", 

(M<f )+(C) = (A? )-«) V^(C), (GS-or ( € E<" (6.14) 

M?(0 = O(j), C^oo, C £ S- or C g 

Given this assessment, undoing the transformation that took us from the RHPs V* - ), (E% K*"), 
(removing the conjugation by m(£) produces MjT* and which solve RHP~^ or RHP'" with new 
asymptotics: 

AC(C) = 0(7)^4= ( \ ! )e-^ CT3 , C-oo, (6.15) 



•C' V2\-i 1 / 

holding uniformly in directions respecting dist(£, E*") > 5 or dist(C,E^) > 5. We show that the 
only conclusion is that M*~(£) = 0. 

6.3.1 Vanishing lemma for RHP^ 

The first step is to fold (and twist) the jumps down to the real line, defining a new matrix Z{£) via 



m 


= Mo 


•(0 




= M - 


'(C) 


m 


= M " 


•(0 


z(C) 


= M - 


•(0 



1 /' 
1 \ / -1 



< arg C < Ivr 



^ 3/2 1 }\ 1 J' |^< ar gC<vr, 



e3 

1 
— es* 5 1 



-7r < argC < — 171", 



■|7r < argC < 0. 



Then Z(Q is an L 2 -solution of the equivalent RHP: 

z+(C) = z.(c)( 4 L ' )■ ^i-x.D], 



h- > 


4,3/2 v 


4 >3/2 


J 


e 3"> 




1 




-1 N 


)• 


1 





z + (() = z_(C)( x I, Ce(0,a], (6-16) 

z+(0 = MO I ° ~ 1 I , CG(«,oo), 
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with now Z(C,) = 0(( 3 / 2 ), compare (|6.15p . Denoting the piece-wise defined jump matrix in (|6.16p 
as V, we notice that, 

0= f Z + {s)Z*_{s)ds= [ Z-(s)V(s)Zt(s)ds. (6.17) 

The first equality holds for any functions in the range of C+ and C_ as may be seen by rational 
approximation. Adding (|6.17p to its conjugate transpose we also find that, 

°=/L z - w l) zl<s)ds+ £ z - (s) { e T l) z "- { ° )d °- 

It is immediate that the first column of Z„ vanishes a.e. on (—00, a], and so {Z\\,Zi\) = 
throughout the lower half plane by analyticity (it lies in the Hardy class H 2 ). From the structure 
of the jump across R one may next conclude that the second column of Z + equals a.e. on (—00, a), 
and by the same reasoning (Z12, Z22) = in the upper half plane. 
If we now set, 

a(C) = Z n (C) in C + and a(() = Z 12 (() in C_ , 
we are led to the the scalar RHP: 

a(C) analytic in C\R_, 

a + (C) = a_(C)ef^ /2 , C € R_, (6-18) 

a(0 = O(r 3/4 ), (eC/l_. 

Note that setting a = Zi\ in C+ and = Z22 in C_ produces the identical RHP . If we can conclude 
that a(() = 0, then Z(() and so Mq(() also vanish identically, proving that ker(I — Cy<-) = 0. 

Lemma 6.2. The unique solution to 16.18]) is a(£) = 0. 

Proof. Up to this point the jump condition has been understood in the sense of 1? . To go further 
it is required that a±(C) are uniformly bounded on the negative real axis. First, since e^ 3 / 2 is 

analytic off R_, one actually has analytic extensions of a+ below R_ and a_ above. That is, (|6.18p 

,3/2 

holds in the sense of continuous boundary values. Further, e^~ decays as C moves below R_ and 

_ r 3/2 ^3/2 

e ^- = decays as as ( moves above R_, and thus both extensions exhibit at least the same 
decay as a(£) itself as (^00. 

Taking the extension of a+(C) into a region ir < arg(£) < tt + e, the Cauchy integral formula 
provides the representation 

«+(C)=/^4^, (6.19) 

in which 

€ = {z : arg(z) = -tt + e/2] U {z : arg(z) = vr - e/2} = <T U £+ , 

3 /2 

oriented counter-clockwise, and a(z) = a(z), a{z)e~ z on (£~, £ + . It follows that for all ( E (—00, <5] 
|a + (£)| is bounded by a constant depending only on 5 > . An identical argument pertains to 

*3/2 

a_(C). To achieve a bound down to ( = 0, we need only note that the jump = e^- for ( < 
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and = 1 for £ > is Holder continuous across zero and the Cauchy transform preserves Holder 
continuity. 

Granted that a(Q is bounded down to R_ from both directions, consider now the effect of 
performing both the above extensions: a+ below to an angle it + ttv/2 and a_ above to an angle 
— 7r — ttv/2 with small v > . The resulting function, denoted by 6(C), can be made to live on a 
subset A of the Riemann surface K formed by gluing together three copies of C cut across R_ in 
the obvious fashion (alleviating the fact that the initial domain swept out a region of angle > 2tt). 

Next bring in the transformation C( w ) = u 2+u which maps the right half of the w-plane onto 
A, taking the positive/negative imaginary axes in co onto the lower/upper boundaries of A. Then 
b(uj) = b(C(to)) is analytic in the open half-plane {uj : 3? (a;) > 0}, bounded in the closed half-plane 
{uj : 3t(u) > 0}, and along the boundary satisfies, < Ce~ c ^ s ^ /2<2+ ^ < C e _c ' s L A theorem of 

Carlson ([19J p. 236) then implies that b = in the right half plane, which is to say that a = 0. □ 

6.3.2 Vanishing lemma for RHP^ 

The steps for RHP^ are mimicked to the point that the needed conclusion hinges on the following. 
Lemma 6.3. The unique solution of the scalar RHP, 

a(C,) analytic in C\R_, 

a + (C) = a-(C)e^- /2+2 ^- /2 , f € R-, (6-20) 

a(C)=0(C 3/4 ), (eC/l_, 

is a = 0. 

Proof. The analysis of (|6.20p is really no different than (|6,18p . The jump is again analytic in C/R_, 
implying that a is continuous and uniformly bounded down to M_ by the same type of extension 
argument. Since a > 0, 

| e | C 3/2_ 2 | 0; | C l/2 | < 1 for ^ . g (7T — £, 7r) U (-7T, -VT + £)} 

and 

|e | C 3/ 2 -2|a| C V 2| < e - cs 3/ 2 for c > Q ^ c = se i( ±7rT£)) s ^ QO) 

thus extensions share the same qualitative features as above. The analogs of b and b are then 
constructed as before and subject to the same conclusions. □ 

7 Properties of the solutions 

We prove a continuity result for M^(£) and M^(£) in the parameter a and establish asymptotics 
of those matrix functions for ± oo; these will lead to Theorems 11.31 and 11.41 
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7.1 Continuity 

The continuity result is based on verifying the condition of the following general fact; see Corollary 
7.103 of [8J for a proof. 

Proposition 7.1. Consider a family of (uniquely solvable) RHP's on a fixed contour, (Yi,v n ), 
n = 1,2,.... Assume the existence of a Voo, such that the RHP (£, v^) posseses a unique solution 
and 

| K - Voo\ |L°°(£)nL 2 (S) ->0, n->CO. (7.1) 

Then, 

||(m n )± - (m oo ) ± || i 2 (S )->0, and \\(m n )(z) - (m oc )(z)| | iO o ( ^)^0, (7.2) 
for n— >oo and any set A which is a positive distance from X. 

The condition (|7.1|) implies that fj, Vn = (I — C Vn )~ l I satisfies \ \n Vn — M^oJ Il 2 (£)^°°- From the 
expression (|6.3p . the statements of (|7.2p easily follow; the first because C± map I? to I? . Further, 



one sees that an estimate of the second type holds for the derivatives, ^[(w-nX-z) — (w-ooX-z)]- This 
is the fact referred to in the proof of Theorem 11.11 

Lemma 7.2. The condition \7.1\ j is satisfied by the RHP's (V^, S*~), the parameter a playing the 
role of n in the Proposition. Continuity holds in each problem down to a = 0. 

Remark 7.3. Note that the problem RHP^ has a-dependence in the contour itself via the segment 
[0, a). In this case, for the I? continuity of the boundary values (M~*)± = (M^)±, we show 
continuity in L 2 (S^) where in the segment [0, a) is extended to [0, a') for any a' > a. 

Proof of Lemma \7.S\ for RHP^ . To employ the conditions of Proposition 17.11 a preliminary con- 
jugation is made to move the dependence of the problem on a > from the contour into the 
jump. 

Consider first the continuity at a point a > 0. In this case the conjugation is affected by the 
same parametrices used in the poof of existence. Set aside a neighborhood of a, U a>E = (a — e, a + e) 
for e > with e <C a. Fix also positive a and b with a < b < a — e and disks D a , Df, enclosing 
a + e (Df, C D a ). Within D a , and for any /3 S U a ,e, we have the parametrices, 

( I JL/*^*' 

Pp(Q = 1 2ni J a *-C aS 
\ 1 

That is, Pp(C) satisfies the jump condition across a < ( < (3. Next define, 

M, {t ,)4 M ^ )m) - ieDb ' (7.3) 
\M~(Q, (€C\D,,. 

Now for all j3 in the defined range we have a family of RHP's on the same contour, with the 
dependence of (3 occurring only in the jump 




30 



Also, for all £ S dDb except ( = b, 

\c p (0-c a (o\ < 



-4/35- 



3/2 



/3 



C 



-(is 



a 



(7.4) 



there being a positive distance £ separating and the interval between a and (3 in this case. On the 
other hand (Cp)±(b) = (C a )±(b). Thus, the above L°° estimate leads to an L 2 (dDb) of the same 
order. 

It follows see that Ma satisfies the criteria of Proposition 17.11 and so M~*(() is continuous at 
any a > in the sense of its boundary values in L 2 (YT' D (C\D(,), and also in L°° for £ exterior 
to Db and away from E~\ This already gives the type of continuity claimed in Theorem 11.31 and 
Claim 15.21 used in the proof of Theorem 11.11 

We complete the analysis by showing the boundary data is L 2 -continuous in the interior Db- 
Inverting the move (|7.3p we have. 



ikT(C) 



M p {C,)Pp\Q 
-Mp{QC p {Q 



(7.5) 




Consider the i-limits of the right hand: we want to show they are continuous in L 2 [b,a + e] as 
j3 ranges in U a ^ £ . First, Mg(£) is analytic inside of with continuous boundary values along 
dDb, excepting the point £ = b. It therefore lies in L°° of that interval, and the conclusions above 
include that j3— >Mp{C,) is continuous in L 2 [b, a + e]. A look at the second line of (|7.5p explains that 
it remains to show that (Cg)± are continuous in L 2 [b, a + e]. But, taking f3 { a from above without 
any loss of generality, 



\\{Cp)±-{C, 



a)±\\L 2 {b,a+e] 



e 3 s ds + 



4vr 2 



ra+e 


r _A s 3/2 

f p e 3 s 




/ ds 


h 


J a S-t 



dt->0. 



Here, the Holder continuity of e _4//3s3/2 produces the vanishing of the integral over [b, a], and the 
first term on the right follows from ||C± o \ \l 2 ^ II ° Hl 2 - This completes the proof for a > 0. 

Turning to the case a = 0, the first step is to delay the jumps to the left of the origin in the 
original problem (V~*, E"*) by considering the equivalent RHP: 



Ce [0,/?], 
Ce [-1,0], 
Ce (-oo,- 




CG{C:arg(l+C)=±|vr}, 
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where M(Q is otherwise analytic and equals (I + 0(( 1 ))m(C) as £— >oo (recall (|6.4p ). This problem 
is obtained from (y~*,T,~*) by setting 

M(C) = M-(0 ( ! 3/2 ° V (ew±, 

y ±es^ 1 J 

where W± is the intersection of C+ or C_ with the region bounded between the rays 

{C : arg(l + C) = ±^vr} and {C : arg(C) = ±^vr}. 

Proving we have I? continuity here for /3 J. will imply the same for the original problem. 
Similar to above, we now set 



where 



W= „ ---- = „ 7' ■ (T.6) 



/(«) = 1 for - 1< s < and /(a) = e"t s3/2 for < s < (3. (7.7) 



Again, the point is that Pp(C) satisfies the jump condition 

(P„)+(0 = (P/3)-(C) ( J /( x ), C € (-1,0). 

Conjugating out by Pp(C) inside a disk -D1/2 = {C : |-z — 1/2| < 1/2} has three affects. First, a 
new jump of P^" 1 ((") is produced along dD^/i- Second, the jump which M(Q has across [0,(3] is 
eliminated. Third, the jump across [—1/2,0] now reads 

/ 1 (CfsUC) \ ( e^T 1 W 1 -(C P ) + (C) \ 



(7.8) 

4,3/2 4,3/2 4 ,3/2 \ 

sC- l + (^)_(C)e3C- -(^) + (C)e-3 C + \ 

4,3/2 

e* c - / 



10 e-!^ /2 1 



e 



We already understand that the jump P^ 1 (C) is continuous in L°° n L 2 of <9Lq/ 2 - To check that 
the jump (|7.8p satisfies the like conditions over —1/2 < £ < 0, note that it is only the (2, l)-entry 
which requires investigation and that term (neglecting the constant 1) may be rewritten as in, 

4,3/2 4,3/2 

(CpU()e^- -(Cp) + {Qe-& (7.9) 

= (CpUO - {c p )+(0 + (CpUQ (^- /2 - 1) - (q,)+(C) (e-t^ /2 - l) 

= -1 + (CpUQ (et^- 2 - l) - (Cp) + (() (e-^T - l) . 



The continuity in L [— 1/2, 0] as (3 [ follows by a computation similar to (|7.6p and the boundedness 

,4,3/2 

of e 3*>± — 1. As for the continuity in L°° recall that the maps C± maintain Holder continuity, 
so there is no problem for £ in the interior of [—1/2,0]. The potential issue of the logarithmic 
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4^3/2 

singularity of (Cp)±(0) as f3 [ is countered by the fact that e^3^=p — 1 vanish to higher order at 
the origin. 

The criteria (|T. 1[) has thus been checked for the new RHP created by the conjugation by Pfl(C) 
defined in (|7.6|) and (|7.7|) within a neighborhood D1/2 of the origin. It remains to invert this move 
and show that L 2 -continuity at a = of M± (and so M±) follows suit. However, the needed 
argument is identical to that given above in (|7.5p and surrounding discussion. □ 

Proof of Lemma for RHP'" . The verification of the conditions in this case is straightforward 
on account of the contour being independent of a from the start. For any positive a and (3, the 



difference of VZ~ and of course vanishes on P_, while on the lines 7± = {C : ar §(0 = if 71 "}) 



Due to the decay of e 4//3< = 3,/2 along 7 =t it is plain that 

|| e3 C - e 2a( > )||L- (7 ±)nL2( 7 ±)^0, 

as /3—+a, the case of a = and j3 [ being no different. □ 

7.2 Asymptotics as ± oo 

As a— > + oo, it is intuitive that the (unique) solution of RHP^ should converge to the solution 
Pa(C) °f t ne RHP defined by the jump conditions, 

(p a ) + (C) = (Pa)-(C) 

(P A )+(C) = (Pa)-(C) | Cm ; ), argC = ±l7r, (7.10) 

(p A ) + (C) = P-(C)( 




with P(C) having the same asymptotics as M - *■(£) as Q— >oo. As is well known, P(C) is given 

2 

explicitly in terms of the Airy function Ai(£) and its derivative . In particular, with uj = es m , let 

Ai(C) Ai(c^ 2 C) ^ 
Ai'(C) AiVO 



no =[ :r 2 s ; i, cec + , 



P(o = I Ai(c) C G 

' Ai'(C) -AiVC) 



(7.11) 



and T(C) = ( ^ j- Then, 



P A (C) = v^e-^PfOed^-f)^, -fvr < argC < fvr, 

P A (C) = v^e— V12 p( C ) e (§C 3/2 -f )oa T(C)- 1 , |tt < arg ( < it, (7.12) 

P A (C) = v / 2^e-^/ 12 P(C)e(l^ 3/2 -f ) CT3 T(C), -ir < argC < -|tt 
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We have the following. 
Lemma 7.4. As a—> + oo, 

M-(C)(Pa)-\0 = (/ + O(e-§ o3/2 )), (7.13) 

uniformly for £ supported away from C\[a, oo). 

If instead — oo, one takes advantage of two facts. First, the jump for RHP^ along arg£ = 
±§7r satisfies 

( e |c 3 / 2 +2«C 1 / 2 1 ) = ( e -2MC 1/2 (i+o(i)) i V for C = o(|a|) ' 
and, second, the unique solution of the RHP: -Pb(C) analytic in C\S~ > , 

(Pb)+(C) = (ft)-(C) ( JL ° )> CG7 ± 



e- 2 ^ 1/2 1, 

(ft)+(C) = (ft)-(C) I ^ J J, CG(-oo,0), 



(7.14) 



Pb(C) = C -3/4 1 f 1 » U/ + 0(|C|- 1 / 2 )V C-oo, (7.15) 



and 

y j l 

is known explicitly in terms of Hankel functions. Here 7 =t are any rays (eventually straight) rays 
extending above and below the negative real axis as in Figure [6l We have in fact already seen the 
solution in part. Set Q(£) to be as defined in (|6.1ip . (|6.12p and (|6.13|) but in regions /, II, and 
III respectively (see again Figure [6|). Then, 

P B (C) = v^Q(C)e- fl/2<T3 . (7.16) 

That (|7.16p satisfies the jumps (|7.14p is immediate from the jump relations for Q(Q- Note that 
replacing the straight lines arg£ = ±27r/3 with (67^ has no affect: Q(C) is analytic off M~ and 
the jump contours may be deformed to accomodate this change. Lastly, the asymptotics (|7.15|) can 
be verified from substituting the formulas, 

h£\0 = y^-i )(i + o(c 1/2 )), H^io = ^-^-f )(i + o(c 1/2 )), 

for £^00 ([l], formulas 9.7.1 - 9.7.4) into the definition of Q(£). 
The analogue of Lemma 17.41 can now be stated. 

Lemma 7.5. Set 

£a(C)=(M-oC) e"- CT3 , for ICI <e\a\, (7.17) 



3 

with any e < 1. Then 7 =t in ^7.14\ ) ma V be chosen in such a way that 

2 



AT-(C) 
uniformly on |£| < e\a\. 



E a (C) P B { C(\a\ - ^0 2 



I + OQal" 1 )), a^-00, (7.18) 
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Figure 6: The contour for Pq and corresponding transformation. 



From Lemmas 17.41 and 17.51 Corollary 11.41 is read off immediately from the explicit forms of Pa 
or E a PQ- we have for instance, 

(M-)u(C) = v^Fe-^AiCOef^ 3 (1 + 0(e~h 3/2 )) 

and 

(M-)n(C) = V^e— /\\a\ - ^C) 1/2 h(Vl{W\ ~ lo)^^ U + OCM -1 )). 

Lemmas 17.41 and 17.51 themselves follow directly from checking condition (17.191) of the below propo- 
sition, the proof of which may be found in [12], Section 7. 



Proposition 7.6. If for a family of L 2 -solvable RHP's (E, v n ) there is the estimate 

C . . 

IK - -f||L°°(E)nL 2 (£) < — , (7-19) 

for a fixed constant C and all large n, then ||C'u„||l 2 (s;)^.I/ 2 (e) = ^(n) ; an< ^ ^ e so ^ u ti° ns satisfy 
m n(C) = I + ^(n) un if orm ly f or C a positive distance from S. 

The fact that (|7.19|) implies a like bound on the operator norm of C Vn actually implies the 
existence of a full asymptotic expansion of M^(() and M*~(£) in powers of e~° 3/2 or a" 1 with 
sectionally analytic coefficients. This is not pursued here. 



Proof of Lemma \7.4\ Define 

p(c) = M-(c)(p A r 1 (c), 

which solves the RHP: 

R(C) analytic in C\[a,oo) 

R + (C) = R-(C) (Va)-(C) ( J ~ 6 'f J (^a): 1 ^)^) , C G [a, oo) (7.20) 
i?(0 = / + O(i), C^oo. 
The jump matrix along [a, oo) can be simplified as in 

(PA) - (C) 1 1 J (PA) - (0 - H (Ai'(C)) 2 Ai(C)Ai'( C ) J 

= i-v r (0- 
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Next, noting the asymptotics, 



Ai(C)= %^e-K 3/2 (i + o(icr 3 / 2 )), 



2 



(see pQ, p. 446), we have that both HVrIIloo^oo) and 1 1 Vr\ \i 2 [a,oo) are bounded by constant multiples 
of e~° S/2 , and the claim follows. □ 



Proof of Lemma 7.5. First we consider the scaled RHP for 

M w (w) = Af*"(|a|w), 

which has the jump conditions: 

( 1 \ 2 / U I 



e 



|a| 3 / 2 (2«; 1 /2_| w 3/2 ) 1 ,,^ 6 --^ 3 -., — 1 _^ ^ 



argu; = ±— 7r, and , u> G 



We will now extract a local parametrix in a neighborhood of w = 0. For w 6 f/ e = {|u>| < e} and 
e < 1 define 

2 

7/ = r](w) = w(l w) 2 

3 

Clearly, r](w) takes C/ £ in a one-to-one fashion onto a open neighborhood C/ £ of v = 0, sending the 
negative real line to itself and the segments argu> = ±2-7r/3 onto rays 7 =t C U £ lying above and 
below the real axis. Extending 7 =t to 00 (smoothy) along straight lines outside of U £ we have the 
jump relations 

( l-2\a\'W l ) ' V E 7± ^ ( ) ' 77 G ^ 

This identifies the choice of j^, and with this choice the above problem is solved by Pedal 3 ??). 
What is the same, Pb(|ck| 3 ^(^)) satisfies the jump relations for M^\w) restricted to \w\ < e (in 
which the upper and lower contours are pulled back to avgw = ±2ir/3. 
Next we perform a second transformation, setting 

f M ( - 1 \w)(m(\a\w)) , \w\ > e, 

mW{w) = { V '} '> ! . (7.21) 

( MW(w)( y E a (\a\w)P B (\a\ i 7](w))j , \w\ < £ 

The definition of E a (\a\w) is given in (|7,17p . It is analytic in |^| < e and so E a (\a\w)PB(\a\ 3 r](w)) 
also shares jump conditions with M^(w) in \w\ < e. 

The point is that M^(w) satisfies a new RHP with jump contour consisting of three pieces: 
the rays = argu; = ±27r/3 n{|u>| > e} and the dU £ the boundary of the disk of radius e. 

On either of the first set of contours, T^, the jump matrix is bounded as in 



m(l«M ( e _ H 3/ 2( l_2 w 3/ 2) I (m(\a\w)) 1 



</ + e -§H 3 H 3/2 



a 2 \z\ x l 2 1 
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which is to say it is I + 0( e -°'«H 3 ) in L 2 n L 00 ^). On dU £ the jump matrix is E a (\a\w) 
PB(\a\ 3 r](w)) (m(|a|w;)) _1 and we compute: for = e and \a\ 3> 1, 

(\ (73/4 

= l + 0{a- 1 ). 

It follows that 

M^ 2 \w) = I + ©(a™ 1 ), uniformly for u; supported away from dU £ UT + UT. 
Undoing the transformations inside of U £ establishes the claim for M^(Q. □ 
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